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Specific heat measurements constitute one of the most powerful experimental methods to probe 
fundamental excitations in solids. After the proposition of Einstein’s model, more than one century 
ago (Annalen der Physik 22 , 180 (1907)), several theoretical models have been proposed to describe 
experimental results. Here we report on a detailed analysis of the two-peak specific heat anomalies 
observed in several materials. Employing a simple multilevel model, varying the spacing between 
the energy levels = {Ei — Eq) and the degeneracy of each energy level Qi, we derive the required 
conditions for the appearance of such anomalies. Our findings indicate that a ratio of A 2 /A 1 ~ 
10 between the energy levels and a high degeneracy of one of the energy levels define the two- 
peaks regime in the specific heat. Our approach accurately matches recent experimental results. 
Furthermore, using a mean-field approach we calculate the specific heat of a degenerate Schottky-like 
system undergoing a ferromagnetic (EM) phase transition. Our results reveal that as the degeneracy 
is increased the Schottky maximum in the specific heat becomes narrow while the peak associated 
with the EM transition remains unaffected. 

PACS numbers: 64.60.A-, 65.40.Ba, 65.60.-l-a, 71.27.-l-a 


I. INTRODUCTION 

In the field of condensed matter Physics, specific heat 
measurements can be considered as a pivotal experi¬ 
mental technique for characterizing the fundamental ex¬ 
citations involved in a certain ph ase transition. In¬ 
deed, phase transitions involving spirP^, charg^, lattic^ 
(phonons) and orbital degrees of freedom, the inter¬ 
play between ferromagnetism and superconductivitj®, 
Schottky-like anomalies in doped compound^, electronic 
levels in finite correlated system^, among other features, 
can be captured by means of high-resolution calorime¬ 
try. Furthermore, the entropy change associated with a 
first-order phase transition, no matter its nature, can be 
directly obtained upon integrating the specific heat over 
T, i.e. C(T)/T, in the temperature range of interest. In 
his seminal paper of IQOiP, in order to explain the devia¬ 
tion of the specific heat of certain materials, like silicon, 
boron and carbon, from the Dulong-Petit’s law, Einstein 
proposed a model based on the assumption that all atoms 
in a solid vibrate independently from each other with the 
same eigen-frequency. The well-known Einstein’s expres¬ 
sion for the phononic specific heat at constant volume, 
Cph,v{T), is given by: 

pOhv 

C^uAT) = ( 1 ) 

hereafter all extensive quantities are defined as per parti¬ 
cle, except when otherwise indicated; f3 = I/T, (fc^ = I 
hereafter), tz is the eigen-frequency of the oscillator and 
h the Planck’s constant. The quantity hv is the so-called 
Einstein temperature 0^. It is worth mentioning that 
recently, one of us made use of Einstein’s model to de¬ 
termine the eigen-energy of the counter-anions libration 
modes in a molecular conductoi!^. However, thought the 
success of the model proposed by Einstein, the theoretical 
description introduced by P. Debye in 191^^21^ constitutes 
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FIG. 1. Main panel: Specific heat (C) as a function of tem¬ 
perature for an hypothetical Debye’s solid (cf. Eq.[^ with 0 _d 
= 120 K. Inset: C/T versus T showing a maximum centered 
at T ~ O.28-0D. 

the hallmark in the description of the phononic contribu¬ 
tion to the specific heat in solids, see e.g.^^. Essentially, 
Debye’s model is based on the hypothesis of a continuous 
isotropic solid. The dispersion relation is linear, i.e. the 
sound velocity is constant and isotropic, being a set of 
eigen-frequencies allowed to the oscillators. In the frame 
of Debye’s model, the phononic specific heat at constant 
volume, ^^ ^(T), reads: 

where x = Qd/T and 0 d is the so-called Debye tempera¬ 
ture. The behavior of the specific heat as a function of T, 
according to Eq.[^ well-known from textbooks, is shown 
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in the main panel of Fig.[^ A particular behavior, usu¬ 
ally not explicitly shown and discussed in textbookJi^, 
is observed by plotting CjT versus T, cf. inset of Fig. ^ 
The maximum centered at 0.28 -©d corresponds to the 
inflexion point of Eq.[^ and manifests itself as a direct 
effect of the accessible energy levels upon increasing the 
temperature. In other words, as the temperature of a 
certain Debye’s solid is increased the number of accessi¬ 
ble energy levels is reduced and, as a consequence, the 
entropy variation rate is reduced. The latter indicates 
that for any system which obeys a Debye-like behavior, 
0_D can be directly estimated by plotting C/T versus T. 
Similarly in Einstein’s model (Eq.[^, C/T versus T has 
a maximum centered at O.38-0£:. A combined system of 
Debye and Einstein phonons can be studied provided the 
Einstein and Debye temperature scales are not so clos^i^. 

At low temperatures, namely for ordinary metals T ~ 
0£)/5O ~ liquid ^He temperature, the following relation 
is valicP^l 

C,{T) =-fT + AT\ (3) 

where 7 is the Sommerfeld coefficient and A = 
127r^/50|). Here the Debye temperature can be di¬ 
rectly obtained via A parameter, whereas the effective 
mass m* of the charge carriers in a metal can be esti¬ 
mated by means of the Sommerfeld coefficienlPi^. The 
latter can be considered the smoking gun, for instance, 
when discovering materials with heavy-fermion-like be¬ 
havior. Moreover, still considering the relevance of high- 
resolution calorimetry, combining Cp{T) and the linear 
thermal expansion coefhcientlSl, ai{T), and making use of 
the Ehrenfest theorem, see e.g.^S^, the uniaxial-pressure 
dependence of the critical temperature for pressure ap¬ 
plied along the Taxis for a second-order phase transition 
can be directly obtained, as follows: 



FIG. 2. Schematic representation of the energy levels Eq, Ei 
and E 2 considered in the model (Eq.[^. The excitation gaps 
Ai and A 2 are dehned as the energy difference relative to the 
ground state energy (Eq). See discussion in the main text for 
details. 

In such a system, the so-called Schottky anomaly man¬ 
ifests itself as a shallow maximum in the specific heat 
data as a function of temperature. 

After this brief introduction, recalling some funda¬ 
mental aspects related to specific heat measurements in 
solids, well-known from text books, see e.g.^^^, below we 
present a multilevel model to describe systems which 
present multiple peaks in the specific heat data. Interest¬ 
ingly enough, such a simple model is capable to describe 
experimental electronic specific heat results of correlated 
electrons systems. 


dT, 

dP, 


= V„ 


,ol 


■Tr_- 




( 4 ) 


II. MATHEMATICAL MODELING 


where Vmoi is the molar volume, and AC refer to 
the thermal expansion and specific heat (per mol) jumps 
at the transition temperature, respectively. The index i 
refers to the crystallographic direction, along which pres¬ 
sure is applied. Rigorously speaking, the Ehrenfest rela¬ 
tion is applicable only for mean-field-like phase transi¬ 
tion, where both Aa and AC present step-like behavior. 
Note that the Ehrenfest relation enables us to determine 
the pressure dependence of T^, purely via thermodynamic 
quantities, i.e., dT^/dPi can be estimated without car¬ 
rying out any experiment under application of external 
pressure. 

Eor a general two-level system, separated by an energy 
gap Ai, the specific heat is described by the following 
expression: 


(l-be-/3Ai)2 ’ 


( 5 ) 


In the frame of the canonical ensemble the sum over 
states, namely the partition function [Z) of a system, is 
given by: 


Z = Tr{e-^^)=^e-^^-, ( 6 ) 

n 


where En refers to the energy level, i.e. the eigen¬ 
value of the system’s Hamiltonian As a matter of 

fact, the partition function encodes all information of the 
physical system, so that once we know the energy eigen¬ 
values of the physical system all thermodynamic observ¬ 
ables can be calculated. This in turn is not the case if 
one has the system Hamiltonian in hancP^, being neces¬ 
sary to make its diagonalization to obtain Z and then 
the physical quantities of interest. Eor the sake of sim¬ 
plicity, in order to describe the two-peak like anomalies 
observed experimentally in several materials (see below), 
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FIG. 3. specific heat (C) as a function of: (a) Temperature (T) and degeneracy of the energy level 1, labelled gi] (b) T and 
degeneracy of the energy level 2, g 2 ', (c) T and energy gap Ai; (d) T and energy gap A 2 . The projection of each curve shown 
in the 3D plots is depicted in the various panels. Details are discussed in the main text. 


we propose a multilevel model approach. This Ansatz 
enables us to obtain directly the thermodynamic quanti¬ 
ties of interest. In other words, within such an approach 
we work in the diagonal representation of the Hamilto¬ 
nian H (cf. Eq.[^. We start by considering a three-level 
modeP^, with Ai — (Ei — Eq) and A 2 = {E 2 — E 0 ), being 
Eo, El and E 2 defined as shown in Fig.[^ the partition 
sum reads 

Z = goe-^^^ + -h (7) 

where gi indicates the degeneracy of each energy level 
(z = 0,1,2). The assumption that the energy gaps A^, 


i.e. energy scales for a generic system, do not depend on 
the temperature is quite realistic. In this regard we refer, 
for instance, to the Schottky modeP^, where the energy 
gap Ai separating the two energy-levels is fixed and thus 
temperature independent. 

Making use of Eq.[^ one can calculate the specific heat, 
employing the following well-known relations 

E = -^{lnZ), (8) 

and 

















4 


C = -I3^ 



( 9 ) 


obtaining 


C = 0 - 




( 10 ) 


where i,j = 0 , 1 , 2 and Ag = 0. Note that considering j 
= 1 in Eq. 10 the specific heat for the Schottky model 
(Eq.[^ is nicely restored. 

For the sake of completeness, we calculated the sum¬ 
mations in Eq.[^ and present below the expression em¬ 
ployed in our analysis: 


[5o + + 52e-^‘^2]2 

gig 2 e~^^^^~*~^^^[Ai(Ai — A 2 ) -I- A 2 (A 2 — Ai)] 

[flO + + g^^-l3A2]2 

( 11 ) 

In the present model, the crucial feature refers to the 
relation between the energy scales of the two excitation 
gaps. The latter together with the degeneracies gi of 
the energy levels, define whether the system presents a 
Schottky-like behavior with a single peak/maximum in 
the specific heat or if it will exhibit two or even three 
(in this case a four energy levels is required) peaks or 
maxima. More specifically, one can say that in such cases 
two and three energy scales govern the Physics of the 
system of interest. 


III. RESULTS AND DISCUSSION 


In what follows, employing Eq. 11 we discuss the re¬ 
quired conditions for the emergence of the two-peak 
anomalies in the specific heat. In Fig.[^ we present 3- 
dimensional (3D) plots of the specific heat as a func¬ 
tion of temperature and degeneracy of the energy levels 
1 (Fig-ia) and 2 (Fig.j^b) and, specific heat as a func¬ 
tion of temperature and the relative value of the energy 
levels 1 (Fig.j^c) and 2 (Fig.j^d) as well. The optimized 
parameters used in such set of fits, appropriate to give 
rise to the emergence of a double peak in the specific 
heat are Ag = 0, Ai = 205J/mol, A 2 = 2327J/mol, 

= 1, (71 = 2, g 2 = 4, i.e. in each 3D plot one of these pa¬ 
rameters was varied while the remaining ones were kept 
constant. For instance, in Fig.j^a) T and gi were used 
as the variable parameters. Note that upon increasing 
the degeneracy gi the double peak in the specific heat 
gradually vanishes. A distinct situation is depicted in 
Fig.§b), where the double peak in the specific heat ap¬ 
pears as the degeneracy g 2 is increased. In Fig. I^c) and 
d) we vary the temperature and the energy levels Ai and 


A 2 , respectively. The present findings indicate two dis¬ 
tinct ways of catching a double-peak in the specific heat, 
i.e. either by varying the degeneracy or by increasing the 
energy levels Ai or A 2 . More specifically, Figs.^a) and 
b) indicate that at least one of the energy levels should 
be highly degenerated to give way to the double-peak in 
the specific heat. Similarly to the Schottky-like anomaly, 
this behavior suggest that the high degeneracy contribute 
dramatically to the specific heat over a quite restricted 
T-window. Such a behavior can be easily understood in 
terms of the system entropy: once the number of acces¬ 
sible states is increased and this is the physical situation 
of a highly degenerated generic system or energy level, 
the entropy of the system is increased and, as a conse¬ 
quence, a double-peak structure shows up in the specific 
heat, since the entropy of the system is proportional to 
the area and can be estimated via AS = C/TdT, 
where Ta and indicate two distinct temperatures in 
the range of interest. In other words, the system entropy 
is increased at the expense of the satellite peak in the spe¬ 
cific heat. Nevertheless, as discussed above, the specific 
heat is a bulk property. Thus in order to map the energy 
levels of a generic system and its degeneracy, microscopic 
experimental techniques like electron spin and magnetic 
resonanc^i^, for instance, are required. We stress that 
the two-maxima observed in the specific heat, obtained 
from the above-discussed Ansatz, is a universal signa¬ 
ture of any system with two distinct characteristic energy 
scales that differ from each other at least by roughly one 
order of magnitude. In fact, such features in the spe¬ 
cific heat are related to various ongoing topics of interest 
in condensed matter physics. In particular, 4/-electrons 
based magnetic systemsi^ show spectacular properties in 
this regard. 

Hereafter we focus on examples of various materials 
classes, in which a two-peak like anomaly in the spe¬ 
cific heat versus temperature have been observed exper¬ 
imentally. We start with the heavy-fermion compound 
Ce 3 Pd 2 oSi 6 , a magnetic field-induced quantum critical 
systenJ^. For this material, J. Custers et al. observed 
the presence of two distinct peaks centered at Tq ~ 
0.5 K and Tn ~ 0.25 K under an external magnetic field 
of 0.5 T in the electronic specific heat. These authors 
attributed the features at Tq and Tjv to a antiferro- 
quadrupolar magnetic order of the Ce 4/ orbitals located 
on the 8c site and to an antiferromagnetic order, respec¬ 
tively. Yet, the authors point out that this scenario is 
compatible with the Tg quartet and Tr doublet ground 
states due to the difference of the crystal-field splitting 
of the Ce atoms on the 4a and 8c lattice site^^. The 
small energy difference of only AT ~ 0.25 K ~ 20/ieV 
between the transition temperatures Tq and T^ sug¬ 
gests that for Ce 3 Pd 2 gSi 6 a subtle change in the sys¬ 
tem total energy is enough to change the character of 
the Ce 4/ orbitals. Interestingly enough, the suppression 
of Tq gives rise to a magnetic field induced quantum 
phase transition. Also for CeAuGe, CeCuGe, CeCuSP^, 
TbPOj^, A 3 Gu 3 (P 04)4 (A = Ca, Sr) and Cu(3- 
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ergy scales A^, employed in our model, give us indirect 
hints of the physical phenomena associated with the spe¬ 
cific heat maxima. In general terms, the description of 
the nature of a phase transition, which usually manifests 
itself as sharp maximum in the specific heat data (see 
Fig.|^ in the low-temperature window, as well as the 
physical mechanism responsible for the broad maximum 
in the high-temperature range require the combination of 
complementary experimental techniques. Yet, it is worth 
mentioning that on the lower temperature flank, namely 
around the phase transition critical temperature T^, the 
discrepancy between fits and the experimental data can 
be attributed to the absence of critical fluctuations, see 
e.g. Refs.l^^HaUg^jjj^ references cited therein, not negligible 
upon approaching T^, in our Ansatz. In the next section 
we discuss a degenerate Schottky model. 


FIG. 4. Electronic specific heat experimental data Cif (cir- 
cle s) a nd theoretical fits (solid lines) for CeAuG^^ using 
Eq.|^with ks = 8.314J/(mol.K), go = 1, gi = 2, g 2 = 4, 
Ai = 205 J/mol (282 J/mol) for the solid pink line (solid blue 
line) and A 2 = 2327 J/mol. 


Chloropyridine)2(N3)J^, Cu(en)2Ni(CN)J^, URu2SiJ^ 
where a magnetic transition followed by the appear¬ 
ance of su per conductivity is observed, (Cei_a:Laa;) 3 AP^, 
PrOs 4 Sbi 2 p^and UPtsP^with two superconducting tran¬ 
sitions, double peaks in the specific have been reported. 
We now present the results of fits that can be compared 
to experimental data set for the specific heat of a partic¬ 
ular heavy fermion compound. In Fig|^we show exper¬ 
imental literature results of the electronic specific heat 
for CeAuG^^ together with theoretical fits employing 
Eq.[TTl It is worth mentioning that the ordinary phononic 
contribution to the specific heat were subtracted by the 
authors of Ref.l^, being thus shown in Fig|^the bare elec¬ 
tronic specific heat contribution originated from the 4/- 
electrons (Cif). Interestingly enough, despite the non¬ 
interacting gas picture we have adopted in this work, 
the present approach is capable to describe an elec¬ 
tronic system, in which correlation effects are present as 
discussed above. In order to determine the nature of 
the electronic/magnetic “excitations” responsible for the 
emergence of a double-peak anomaly in the specific heat 
data of CeAuGe, as well for the various materials above- 
mentioned, microscopic/spectroscopic data are required 
and constitute a topic out of the scope of the present 
work. It is worth mentioning that inelastic neutron scat¬ 
tering investigation^^, carried out on CeAuGe at 15 K, 
revealed the presence of pronounced crystal field excita¬ 
tions at 24.3 meV (~ 282 K), which roughly corresponds 
to one of the energy scales, namely Ai = 205 J/mol, 
employed in one of the fits (solid pink line) depicted 
in Fig.|^ Note that a fit using Ai = 282 J/mol (solid 
blue line in Fig.|^ does not deviate so much from the 
experimental data. Hence, we figure out that the en- 


IV. DEGENERATE SCHOTTKY MODEL 

Following the proposal of modelling and understanding 
of the double peak anomalies in specific heat data, in this 
Section we focus on the description of a degenerate Schot¬ 
tky model. To this end, we consider a two level system 
(ei = —|A,e 2 = |A, with respective degeneracies ^ 1 ^ 2 ) 
plus a spin ^ Zeeman splitting contribution i/isiJ, cf. 
scheme shown in the inset of Fig.[^ The partition func¬ 
tion and equation for the energy (obtained using Eq.[^ 
per particle is, disregarding an additive constant, given 
respectively bji^ 




( 12 ) 


E = —mH — 


Ae-^^ 

{1 + ne-P^)' 


(13) 


where H = H 1 /H 2 , H is an external magnetic field and 
P refers to the inverse temperature. The magnetization 
m and the specific heat per particle are in turn given 
respectively by: 


m = fiBtanh{P^BH), (14) 


dE d , niPAfe- 

~ dT~ (l-bHe-/3A)2' 


(15) 


For an ordinary paramagnet, cf. discussion in Section I 
(see Eq.[^, the magnetic contribution to the specific heat 
is Schottky type, namely it is described by the second 
term of Eq.[T^with H = 1 and an energy gap A = 2^bH, 
so that we have two decoupled Schottky like systems. 
We assume now that a ferromagnetic phase transition. 
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T(K) 


FIG. 5. Main panel: specific heat as a function of temper¬ 
ature for various energy gaps A and degeneracy ratio SI = 
cf. indicated in the label. Inset: energy scheme of the 
degenerate Schottky model. See details in the main text. 


units) and degeneracy ratio SI, namely A = 10 and SI = 
5 (black solid line); A = 20 and SI = 5 (blue solid line); 


A = 10 and SI = 10 (red solid line) employing Eq. 17 


The low-T peak, with a mean-field like shape, remains 
unaffected upon tuning A(>> 1) and SI. The Schottky 
maximum, however, becomes broader and shifts to higher 
temperatures as A is increased, being the peak inten¬ 
sity also increased as SI increases cf. expected for a two- 
level model, a feature that can be understood in terms 
of entropy arguments as discussed in a previous section. 
It is worth mentioning that between the ferromagnetic 
transition and Schottky maximum is exponentially small, 
i.e. the ferromagnetic phase transition produces a jump 
in the specific heat data on top of a non-magnetic contri¬ 
bution. For A values close to unity (not shown in Fig.[^ 
the high-T maximum tail evolves into the low-T, being 
the sharp mean-field like jump AC = 3 at T = 1 pre¬ 
served. Thus, our analysis show that a simple Schottky 
model coupled to a ferromagnetic phase transition is also 
amenable to fit emerging double peak/maxima features 
in measured specific heat data. 


which is a typical many body effect, takes place at a 
certain critical temperature Tc within the Curie-Weiss 
molecular mean-field modeP^. We introduce then the 
macroscopic molecular field phenomenological parameter 
A and an effective magnetic field (iJe//), as follows: 

HeS = Hext + (16) 

At zero external magnetic field m is calculated in a 
self-consistent fashion via the following equation: 

TO =/is tanh(/3/iBAm), (17) 

being thus now the magnetic gap A = 2(1^11 = 2fj,BXm 
a function of the order parameter to. Scaling the mag¬ 
netization to /tb, the energies (and temperatures) to Tc, 
with Tc = XfJ-B, we have: 


V. CONCLUSIONS 

To summarize, we have employed a multilevel model to 
describe specific heat data as a function of temperature 
with two maxima. Our results suggest that: i) a ratio 
of A 2 /A 1 w 10 between the energy levels and ii) the en¬ 
ergy levels degeneracy, govern the two-peaks regime in 
the specific heat. Apart from the nature of the entities 
responsible for the emergence of the maxima in the spe¬ 
cific heat data, this simple model describes nicely recent 
literature results for the CeAuGe compound. Further¬ 
more, the specific heat of a degenerate Schottky model 
was calculated. Our analysis demonstrate that double 
peak/maxima features in specific heat data can also stem 
from a two-level model coupled to a ferromagnetic phase 
transition. We trust that our findings pave the way to¬ 
wards the description of experimental specific heat data 
for other systems in both soft and condensed matter 
Physics. 




fI(/3A)^e-^^ 

(l-bfIe-/5A)2’ 


(18) 
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